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f  ABSTRACT 

V 

This  paper  is  a  four  part  dialectic  between  theory  and  experimental  observations  and 
measurements  of  two  dimensional  periodic  surface  waves  in  a  vessel  of  rectangular  cross 
section.  The  first  two  parts  consist  of  minor  contributions  to  the  existing  theory  and 
experimental  measurements  of  elliptic  (deep  water)  standing  waves  and  the  (shallow 
water)  travelling  hydraulic  jump  A  new-  theory  is  proposed  in  the  third  part  which 
explains  the  transition  from  the  elliptic  standing  wave  to  the  hyperbolic  travelling 
hydraulic  jump.  The  result  is  a  family  of  standing  cnoidal  waxes.  To  lowest  order  the 
solution  is  composed  of  two  non-interacting  periodic  progressive  cnoidal  waves  travelling 
in  opposite  directions.  At  higher  order  the  two  wave  trains  interact  but  the  interaction 
produces  nothing  of  qualitative  significance  except  when  the  cnoidal  standing  wave  of 
highest  amplitude  is  considered.  This  family  of  waves  is  validated  by  experimental 
observation  in  part  four.  In  addition  in  part  four  observation  of  compound  cnoidal 
standing  waves  leads  to  a  conjecture  about  a  sequence  of  continuous  wave  forms  in 
parameter  space  leading  from  the  linear  sinusoidal  wave  to  the  travelling  hydraulic 
jump.  f 
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ON  THE  TWO-DIMENSIONAL  PERIODIC  SURFACE  WAVES 
OCCURRING  IN  RECTANGULAR  VESSELS: 
THEORY  VERSUS  EXPERIMENT 


Thomas  J.  Bridges 

1.  Introduction 

This  report  is  a  dialectic  between  theory  and  experiment.  It  is  composed  of  four 
parts;  two  on  theory  and  two  on  the  experiment.  The  subject  is  two-dimensional 
periodic  surface  waves  which  occur  in  a  partially  filled  rectangular  vessel. 

The  theory  will  be  concerned  with  an  idealized  fluid:  i.e.  inviscid  and  irrotational 
and  without  surface  tension  at  the  surface.  The  rectangular  vessel  of  length  2 a  contain¬ 
ing  fluid  of  still  water  level  h,  when  a  forcing  function  is  present,  rotates  harmonically 
about  a  fixed  point.  Such  a  configuration  is  depicted  in  Figure  1.1.  The  parameters 
d  and  -y  fix  the  point  of  rotation.  6  is  the  rotary  displacement  at  frequency  u;  of  the 
vessel.  Under  the  conditions  specified  the  governing  equations  are, 


d2(p  d24> 


dx 2  dy 2 


=  0  for  -  a  <  x  <  a, -h  <  y  <  rj{x,t) 


(1.1) 


d0 

—  =  -JIu  on  x 


(1.2a) 


and  on  y  =  r/(x,t) 


dn  d(j>  dr)  dr)  d<j> 

dt  +  dxdx  ^  dx  dy^X 


d<j>  lr d<j>  d<f> .  1  _ ^ .  o  o\ 

m  +  2 [<5>  +  V  ]  +  n(”^  * '*> +  2n + 1  > 

+  T)  gcosO  4-  lid  cos -y  —  fl2dsin/y|  +  xjgsinfl  -  lid  sin  7  -  fl2dcos'yj  =  0  (1.4) 

where  <p  is  the  velocity  potential  and  rj  is  the  surface  displacement,  and  Q  =  0,  the 
angular  velocity  of  the  vessel.  The  derivation  of  these  equations,  under  the  assumption 
of  a  uniform  vorticity  field,  is  given  in  Appendix  I. 

In  the  absence  of  the  forcing  function  the  two  relevant  dimensionless  parameters 
are  6  =  h/2a  and  t  =  H/h  where  H  is  a  measure  of  the  wave  amplitude  above  the  still 
water  level.  The  introduction  of  the  forcing  function  brings  in  the  parameters  d,  7,  w, 
and  d  where  d  is  the  amplitude  of  the  harmonic  forcing  function  0. 

The  surface  waves  may  be  classified  with  the  parameters  £  and  6.  When  (  <<  62 
(which  includes  the  case  £  — >  0,  i.e.  linear  waves)  the  surface  waves  occurring  are 
elliptic  standing  waves.  When  £  >>  <52  the  travelling  hydraulic  jump  occurs.  A  theory 
for  both  these  regions  has  previously  been  established.  A  theory  for  £  <<  6 2  was 
established  by  Tajdbakhsh  and  Keller)  15]  and  a  theory  for  the  region  e  >>  62  has  been 
established  by  Chester[2,3j  and  Verhagen  and  Van  Wijngaarden[l9|.  These  two  theories 
with  improvements  are  given  in  Section  2.  In  Section  3  experimental  observations 
corresponding  to  the  theory  in  Section  2  are  contained.  This  includes  photos  and  time 
series  for  elliptic  standing  waves  and  the  travelling  hydraulic  jump.  Both  of  these  wave 
types  have  been  studied  and  observed  in  the  past.  However  the  region  in  between  has 
recieved  almost  no  attention. 
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In  Section  4  a  theory  is  proposed  for  surface  waves  in  the  intermediate  region 
when  <  =  0(6  2).  The  result  is  that  a  family  of  cnoidal  standing  waves  exists  in  this 
region.  Then  in  Section  5  experimental  verification  of  these  waves  is  discussed  along 
with  another  class  of  compound  standing  waves.  It  is  conjectured  that  these  compound 
standing  waves  are  combinations  of  various  mode  cnoidal  standing  waves,  and  as  the 
amplitude  is  increased  higher  mode  cnoidal  standing  waves  are  excited  and  eventually 
there  is  enough  harmonics  to  produce  a  travelling  hydraulic  jump.  If,  in  fact,  the  family 
of  cnoidal  standing  waves  are  the  only  eigenfunctions  existing  in  the  intermediate  range 
between  the  elliptic  standing  wave  and  the  hydraulic  jump  the  the  following  scenario 
can  be  constructed  for  the  complete  range  of  t,6  when  6  is  sufficiently  small. 

When  t  <<  S2  (and  6  small,  including  the  limit  as  the  amplitude  goes  to  zero, 
i.e.  linear  waves)  an  elliptic  standing  wave  forms  which  bifurcates  from  the  state  of 
rest.  This  wave  is  well  described  by  the  theory  of  Tajdbakhsh  and  Keller.  It  is  a 
weakly  nonlinear  wave  containing  only  a  few  significant  harmonics.  However  this  wave 
breaks  down  as  <  approaches  62  in  magnitude.  This  may  be  shown  with  the  solution 
of  Tajdbakhsh  and  Keller  itself  but  in  an  interesting  piece  of  work  Vanden-Broeck  and 
Schwartz  18]  using  a  numerical  scheme  have  shown  that  the  bifurcation  of  the  natural 
frequency  has  a  turning  point  when  c  =  0(6 2)  amd  the  theory  of  Tajdbakhsh  and 
Keller  can  not  get  around  the  turning  point.  Vanden-Broeck  and  Schwartz  computed 
numerical  solutions  for  the  region  outside  the  turning  point  and  conjectured  that  there 
might  be  a  new  class  of  standing  waves  in  this  region.  By  devising  a  perturbation 
scheme  for  small  <  and  6  such  that  the  ratio  e  to  6 2  is  constant,  it  is  found  in  Section  4 
that  this,  in  fact,  is  where  the  family  of  cnoidal  standing  waves  lives.  These  waves  are 
much  richer  harmonically  as  they  have  a  complete  Fourier  series. 
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As  the  amplitude  is  increased  further  the  class  of  solutions  observed  are  combina¬ 
tions  of  the  set  of  cnoidal  standing  waves  selected  by  the  forcing  function.  Observance 
of  this  phenomena  in  the  experiment  suggests  that  as  <  is  increased  beyond  6 2  the  nat¬ 
ural  frequencies  of  the  higher  mode  cnoidal  standing  waves  become  integral  multiples 
of  the  fundamental  in  sequence.  First  the  3-mode  is  excited  and  a  compound  1-3  mode 
cnoidal  standing  wave  is  observed.  As  the  amplitude  is  increased  further  relative  to  62 
a  5-mode  is  added  to  the  wave  form  and  this  forms  the  compound  wiggle  wave  discussed 
in  Section  5.  Further  increase  of  the  amplitude  introduces  higher  and  higher  harmonics 
until  eventually  the  wave  field  jumps  to  the  travelling  hydraulic  jump.  Interestingly, 
this  cascading  process,  unlike  turbulence  for  example,  is  not  dynamic.  At  a  fixed  am¬ 
plitude  each  compound  wave  in  the  sequence  is  (apparently)  stable  as  evidenced  by 
the  results  in  Section  5.  Therefore  this  sequence  of  compound  standing  waves  is  a 
sequence  of  equilibrium  states  in  parameter  space  which  eventually  leads  to  that  flood 
of  harmonics;  the  travelling  hydraulic  jump. 

Obviously  much  work  remains  to  substantiate  these  claims  but  it  does  suggest 
a  continuity  in  parameter  space  between  sinusoidal  waves  and  the  hydraulic  jump. 
However  as  evidence  of  the  lack  of  concrete  and  rigorous  theory  in  this  area  there  has 
not  even  been  an  existence  proof  for  elliptic  standing  waves  when  e  <<  6 2  but  c  still 
finite.  So  much  work  remains  to  be  done  on  this  problem. 


2.  Theory  -  I:  Elliptic  Standing  Waves  and  Hydraulic  Jumps 


2.1  Introduction 

There  are  essentially  two  existing  theories  describing  the  nonlinear  waves  in  a 
rectangular  basin.  The  first  is  the  theory  of  Tadjbakhsh  and  Keller  [15}  which  finds  the 
linear  eigenvalues  and  the  standing  wave  solutions  emitted  by  them.  These  waves  will 
be  referred  to  as  elliptic  standing  waves  since  they  are  governed  by  an  elliptic  partial 
differential  equation.  It  will  be  shown  however  that  these  waves  are  valid  for  only  a 
small  range  of  amplitude,  breaking  down  when  t  =  0(£2). 

The  other  theory  is  that  proposed  by  Chester  [2,3}  and  Verhagen  and  Van  Wijn- 
gaarden  [19].  They  consider  the  region  where  t  >>  62  (shallow  water,  finite  amplitude) 
and  find  that  there  are  no  periodic  solutions  in  this  region  because  dispersion  is  neg¬ 
ligible  here.  But  the  introduction  of  a  forcing  function  renders  the  problem  solvable. 
Periodic  solutions  which  are  discontinuous  may  be  found  in  this  case.  The  discontinu¬ 
ous  solutions  represent  a  family  of  travelling  hydraulic  jumps.  It  is  interesting  to  note 
that  the  hydraulic  jump  is  not  an  eigenfunction.  It  is,  sort  of,  a  state  of  degeneracy, 
where  all  the  standing  wave  modes  coalesce  to  form  this  odd  wave. 

These  two  theories,  however,  leave  a  void  in  the  region  e  =  0(62).  A  theory  is 
proposed  in  Section  4  for  this  region.  It  is  found  that  a  family  of  cnoidal  standing 
waves  exists  in  this  region. 

First,  the  theory  of  Tadjbakhsh  and  Keller  will  be  given  with  some  embellishment. 
Their  analysis  was  for  the  fundamental  mode  only.  Here  it  is  extended  to  arbitrary 
mode  number.  The  main  results  only  will  be  given,  the  details  of  the  procedure  can  be 
found  in  their  paper. 
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2.2  Elliptic  Standing  Waves 


The  independent  and  dependent  variables  are  scaled  in  the  following  way:  x  — »  2a, 
y  h,t  —•  l/w,  w  -»  y/gh/2a ,  r?  -♦  <6,  and  <t>  —  i2as/gh,  where  2a  is  the  tank  length, 
At  the  still  water  depth,  and  the  dimensionless  parameters  are  «  =  H/h  and  6  =  Ai/2a, 
where  H  is  a  measure  of  the  maximum  height  of  a  wave  above  the  still  water  level. 
This  scaling  is  appropriate  for  fluid  of  finite  depth  only  and  must  be  recast  for  the  case 
of  infinite  depth.  The  governing  set  of  equations  and  boundary  conditions,  when  the 
forcing  function  is  absent,  are, 


'0  +  0  =  ° 

d<j>  1 

-2-  =  0  at  x  =  ±- 

<?x  2 

dd> 

—  =  0  at  y  =  ~  1 

dy 

dr]  dr]  d<f>  1  d<f>  , 

"*+‘Ka;-iS3r0 


(2.2a) 

(2.26) 

(2.3a) 


+  (M*) 

To  obtain  a  formal  solution  to  this  problem  in  the  limit  as  t  tends  to  zero,  a  power 
series  solution  in  t  is  sought, 


<f>  =  <j>0  +  2  + 


r]  -  r] o  +  +  t  r]2  +  ■ 

w  =  wq  +  fWj  +  c*W2  ■+• 


(2.4a) 

(2.46) 

(2.4c) 


The  free  surface  boundary  conditions  are  expanded  in  a  Taylor  series  about  y  —  0, 
the  expressions  (2.4)  are  then  substituted  into  the  set  (2.1)-(2.3).  Equating  terms 
proportional  to  like  powers  of  t  to  zero  results  in  the  sequence  of  elliptic  boundary 
value  problems, 


>d2<t>j  ^  d24>j 


0  for  r<x<-,-l<y<0 


at  y  =  0 


(2.7a) 


14,0  dt  6 2  dy  Kj 
d<t>j  _ 

w0  +  t?y  =  5y  at  y  =  0 


(2.7k) 


for  j  —  0, 1, . . ..  The  Rj  and  S} ,  which  are  essentially  the  same  as  those  found  in  the 
paper  of  Tajdbakhsh  and  Keller,  are  dependent  on  terms  of  previous  order  only.  For 
the  linear  problem  R0  =  So  =  0  and  it  is  straightforward  to  show  that  the  zeroth  order 


solution  is 


w°  =  y  tanh  am6 

rj o  =  cos  amx  sint 
cosh  am6(l  +  y) 

00  =  W0  - - - 7 -  cos  amx  cos  i 

cosh  omo 


(2.10) 


where  am  =  mn,  m  =  1,2,...,  and  x  =  x  +  |.  The  points  w0  for  each  m  are 
the  bifurcation  points.  The  higher  order  terms  will  then  give  the  solutions  a'.oii®  the 
branches  emitted  by  these  points.  The  approach  is  straightforward  and  the  details  will 
be  omitted. 


The  first  order  solution  is, 


wi  =  0 


T7 1  =  B ii  +  B 12  cos2f  cos2amx 

. 


,  .  cosh  2am6  (1  +  y)  _  . 

<t>\  =  An  t  +  An  sin  2t  +  A\z - — - - — -  cos  2amx  sin  2< 

cosh  2 Ot frtv 


(2.11) 

(2.12) 


(2.13) 


Obviously,  when  t  —  amb  the  tenet  of  the  expansion  is  violated.  Namely,  higher  order 
terms  should  be  of  higher  order.  This  suggests  that  e m(b)  -  0(b).  However,  applying 
the  same  argument  to  rj  results  in 


lim  rj(z,  t)  -  cos  amz  cost 
6—0 
e  fixed 


Ze 

8  <b* 


cos2omz  cos2t 


-  256- <  COS  3am;C  sin  •  (2  2°) 

Contrary  to  the  frequency,  the  wave  height  expansion  breaks  down  when  c  ~  (am6)2. 
At  this  point  there  is  a  spillover  of  the  higher  harmonics  onto  the  zeroth  order  term, 
again  violating  the  tenet  of  the  expansion.  This  results  in  a  lower  estimate  for  the 
domain  of  validity  in  the  amplitude  of  <m(6)  =  0(62).  It  can  therefore  be  concluded 
that  the  classical  standing  wave  solution  is  valid  for  any  depth  such  that  0  ^  c  ^  f rn  (*) 
with  fm(6)  =  0(62).  In  this  range  the  solutions  are  for  elliptic  waves.  The  boundary 
value  problem  is  essentially  elliptic  with  merely  a  time  dependent  domain. 


The  effect  of  amplitude  on  the  frequency  for  waves  of  this  type  can  be  seen  in  the 
following  bifurcation  diagrams.  Figures  2.1,  2.2,  and  2.3  are  bifurcation  diagrams,  for 
the  frequency,  at  b  =  .4,  b  =  .2,  and  b  =  .1  respectively.  The  bifurcation  is  subcritical 
or  supercritical  depending  on  the  value  of  b.  The  switch  occurs  when  u>2  =  0.  Setting 
■jj2  -  0  results  in  the  equation, 


tanh2  amb  =  *  (29  +  4\/74)  =  4  (29  -  4x^74)  4  -  1 


(2.21) 


Solving  this  shows  that  u>2  =  0  when  amb  ~  1.0581.  Therefore  the  first  mode  switches 
direction  when  b  —  .34,  the  second  when  b  ~  .17,  and  the  third  when  <5  —  .11. 


Figure  2.1  Bifurcation  diagram  for 
the  first  five  two  dimensional  elliptic 
standing  wave  modes  with  6  =  A. 


Figure  2.2  Bifurcation  diagram  for 
the  first  five  two  dimensional  elliptic 
standing  wave  modes  with  6  — 


Figure  2.3  Bifurcation  diagram  for 
the  first  five  two  dimensional  elliptic 
standing  wave  modes  with  6  = 


Figures  2.4,  2.5,  and  2.6  give  examples  of  the  wave  height  for  modes  1,  2,  and  3 
respectively.  They  are  for  6  =  0.4.  The  solutions  are  similar  at  other  values  of  6  and 
£  when  £  <  (m.  The  first  mode  is  the  classical  standing  wave  often  referred  to  as  the 
slosh  mode,  and  the  second  mode  is  the  mode  studied  by  Taylor[l6]  in  his  study  of  the 
highest  standing  wave.  In  Figure  2.7  the  wave  height  at  the  left  tank  wall  is  plotted 
as  a  function  of  time.  This  time  series  is  typical  for  the  elliptic  standing  wave  when 
£  <<  62. 

As  pointed  out  by  Tadjbakhsh  and  Keller  and  Concus  [6],  the  above  expansion 
breaks  down,  even  when  e  <  em,  if  a  higher  mode  becomes  an  integral  multiple  of  the 
fundamental.  This  happens  when 

m  tanh  mnS  —  p 2  tanh  ir6  (2.22) 


for  any  positive  integer  p  and  m  =  2,3, _ When  the  condition  (2.22)  is  met,  an  inter¬ 

nal  resonance  occurs.  In  other  words  higher  harmonics  of  the  fundamental  generated 

u,.  _ i: _ _ ~  _ t: _ 3 _  un _ _ _ _ _ _ _ . _ _ _ : _ 


Figure  2.4  Spatial  distribution  of  the 
wave  height,  for  the  first  mode,  at 
fixed  time  w’ith  6  =  and  e  =  .15. 


Figure  2.5  Spatial  distribution  of  the 
wave  height,  for  the  second  mode,  at 
fixed  time  with  6  =  ^5,  and  t  =  .12. 


Figure  2.6  Spatial  distribution  of  the 
wave  height,  for  the  third  mode,  at 
fixed  time  with  6  - ■A.  and  <  -  10. 
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Figure  2.7  The  time  series  for  the  wave  height  at  the  left  tank  wall  for  the 
first  mode  with  t  —  .07  and  6  =  |. 
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Figure  2.8  The  time  series  of  the  wave  height  for  Chester's  discontinuous 
solution  of  the  forced  shallow  water  equations. 
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Although  these  singularities  appear  to  be  ominous  they  too  are  of  higher  order.  The 
first  singularity,  a 3  =  2o\  will  not  arise  in  the  expansion  until  fifth  order,  at  least, 
which  renders  it  inconsequential.  Also  the  other  two  singularities  above  will  arise  at 
the  fifth  or  higher  order.  Presumably  a  proper  modification  of  the  expansion  could 
be  constructed  (Mack  [lOj  has  done  this  for  the  critical  depths  arising  in  the  circular 
standing  wave  problem,  and  Vanden-Broeck  [17]  has  done  this  for  the  critical  points 
arising  in  the  capillary  standing  wave  problem),  but  it  is  unlikely  that  the  contribution 
to  the  overall  solution  in  this  case  will  be  noticeable. 

In  light  of  these  facts,  the  expansion  of  Tadjbakhsh  and  Keller  will  be  valid,  to  the 
order  given,  for  all  depths,  but  only  for  a  limited  range  of  amplitudes,  e  <  em. 

It  is  worth  noting  that  the  problem  for  elliptic  standing  waves  in  the  presence  of  a 
forcing  function,  where  the  forcing  function  is  near  the  fundamental  natural  frequency, 
has  been  studied  by  Moisiev  [12],  Faltinsen[7],  and  Lou.  et  al.  [9].  The  solutions  are 
obviously  similar  to  the  eigenfunctions  discussed  here,  but  the  important  result  is  the 
relationship  between  the  amplitude  of  the  forcing  function  and  the  response  amplitude. 
Moisiev  showed  that,  when  the  forcing  frequency  is  near  the  linear  resonant  frequency, 
the  response  amplitude  e  =  0(0 */3)  where  0  is  the  amplitude  of  the  forcing  function. 
This  means  that  the  nature  of  the  nonlinearity  of  the  governing  partial  differential 
equation  and  boundary  conditions  (the  operator)  is  cu6ic  in  the  neighborhood  of  the 
bifurcation  points.  It  will  become  clear  in  the  analysis  in  Section  4  that  the  nature  of 
the  nonlinearity  becomes  quadratic  when  e  =  0(62),  and  when  e  >>  62  the  nonlinear¬ 
ity  remains  quadratic.  Therefore  this  operator  has  the  interesting  property  that  the 
nature  of  the  nonlinearity  changes  as  the  amplitude  of  the  response  increases.  These 
conjectures  are  based  on  the  interesting,  and  useful,  theory  of  imperfect  bifurcation  [13], 
and  singular  perturbation  of  bifurcations  [11],  proposed  by  B.  Matkowsky  and  E.  Reiss. 
In  that  work,  they  propose  a  method  for  determining  the  degree  of  the  nonlinearity 
of  operators  in  the  neighborhood  of  singularities  (such  as  bifurcation  points).  This 
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information  is  then  used  to  obtain  the  proper  expansion  parameter  for  the  analysis  of 
the  imperfect  bifurcation  problem  (a  forcing  function  qualifies  as  an  imperfection). 

2.3  The  Travelling  Hydraulic  Jump 

The  other  theory  proposed  for  waves  in  rectangular  basins  is  Chesters  [2,3]  “shal¬ 
low  water  theory”.  This  theory  is  derived  by  taking  the  limit  as  6  — »  0  with  fixed, 
although  small,  amplitude.  To  proceed  the  equations  (1.1)-(1.4)  are  scaled  in  the  man¬ 
ner  described  earlier  except  4>  and  r;  are  now  independent  of  t.  The  forcing  terms  are 
retained  and  their  role  in  the  existence  of  periodic  soluitons  in  this  region  of  the  6  —  c 
plane  will  be  elaborated. 

In  the  limit  as  6  — *  0  the  Laplace  equation  in  conjunction  with  the  bottom  bound¬ 
ary  condition  gives 

4>  =  <t> o  +  n xy62  -  (2.23) 

where  <j>0  is  a  function  of  i  and  t  only.  Substituting  this  expression  into  the  free  surface 
boundary  conditions  and  neglecting  terms  of  0(62)  and  higher  results  in 


d<t>odri  ,,  ^d24>  o  „ 

axax  +  (l+,)^  =0 

(2.24) 

1  (d<j>  0\2  . 

+  2(aJ  +  *>  +  0x  =  ° 

(2.25) 

But  the  substitution  u  =  d4>o/dx  results  in  the  common  form  of  the  shallow  water 
equations 

+uS+(i+,?)£=o  (2-26) 


with  boundary  conditions  u(±  %,t)  =  0.  This  formulation,  and  the  subsequent  analysis, 
differs  somewhat  from  that  of  Chester  but  the  essentials  are  the  same.  The  forcing 
function  is  taken  to  have  the  form  6  =  c2cosf.  The  amplitude  of  the  forcing  function 
is  taken  proportional  to  the  response  amplitude  squared.  Solutions  are  sought  in  the 
neighborhood  of  resonance  and  the  nonlinearity  and  forcing  function  must  balance  in 
order  to  obtain  a  periodic  solution.  This  suggests  that  the  response  in  the  neighborhood 
of  resonance  is  of  Od^l'^2).  It  also  suggests  that  the  nonlinearity  of  the  operator  is 
quadratic  in  the  neighborhood  of  the  linear  eigenvalues.  This  is  to  be  contrasted  with 
the  result  when  the  operator  was  elliptic.  The  governing  equations  are  now  hyperbolic 
as  well.  Therefore,  as  the  amplitude  of  the  wave  increases  the  governing  equations  and 
boundary  conditions  change  type  from  elliptic  to  hyperbolic. 

To  obtain  a  formal  solution  of  this  problem  in  the  limit  as  t  tends  to  zero,  a  power 
series  solution  in  e  is  sought 


W  =  Wo  "f"  +  •  • 

u  =  eui  +  c2U2  + 


*1  —  trll  +  t  1 2  + 


(2.28a) 
(2.28  b) 
(2.28  c) 


Substituting  and  equating  terms  proportional  to  like  powers  of  £  to  zero  results  in  a 
sequence  of  inhomogeneous  wave  equations.  Combining  the  two  first  order  equations 
results  in  a  homogeneous  wave  equation  for  U|, 


U  uj  =  0 


(2.29) 


where  u  is  the  D’Alembertian, 


2  d2  d2 

u~w°dt2  dx1 


(2.30) 


with  homogeneous  boundary  conditions.  Retaining  the  mode  number  as  arbitrary,  the 
general  solution  is 

«,(*,<)  =  /(£)  -  fix)  (2-31) 


where  £  =  t  +  amx  and  x  =  t  -  Qmz,  and  wo  =  am  -  mx.  The  wave  height  haa  a 
solution, 

,i(x,l)  =  -/(0-/(x)  (2-32) 

The  function  /  is  unknown  at  this  order.  It  will  be  found  by  imposing  the  solvability 
condition  at  the  next  order.  The  second  order  equations  are 


dv 2  ,  du2  ,  drn  dr) x  dm 

of  dx  at  ox  ox 

du2  ,  dr) 2  «9ui  ,  A 

w0-r-  +  —  +w  i—  +  uj— —  +  cosf  =  0 
of  ox  of  ox 


(2.33) 

(2.34) 


which  can  be  combined  to  form  an  inhomogeneous  wave  equation  for  U2.  After  substi¬ 
tuting  the  first  order  result  into  the  right  hand  side,  the  result  is 


Uu,  =  -  2amUl|/"(0  -  /"(*)]  -  3«^[(/({)/'(f))'  -  (/(x)/'(x))1 


(2.35) 


There  is  a  solution  to  (2.35)  if  and  only  if  the  right  hand  side  satisfies  a  Fredholm 
alternative.  The  necessary  Fredholm  alternative  for  the  inhomogeneous  wave  equation 
is  given  in  Appendix  III.  Application  of  this  results  in  an  ordinary  differential  equation 
for  /, 

4«^wI/"(x)  +  6<*m(/(x)/'(x))'  +  4amcos(|)  =  0  (2.36) 


which  may  be  integrated  twice  to  yield, 


/2(*)  + £'<*)- cos,f)-c’=° 


(2.37) 


where  e2  =  /Q2am  /2(x)dx.  It  is  interesting  to  note  that  if  the  forcing  function  were 
absent  in  (2.37)  the  solution  would  be  f(\)  —  constant,  but  the  constant  would  have 
to  be  zero  to  satisfy  mass  conservation.  This  shows  that  there  are  no  periodic  solutions 
bifurcating  from  the  state  of  rest  for  the  shallow  water  equations  (2.26)  and  (2.27)  when 
the  forcing  function  is  absent.  In  fact,  Keller  and  Ting  [8]  have  rigorously  proved  that 
for  any  set  of  quasilinear  first  order  set  of  hyperbolic  equations  with  two  unknowns, 
the  solutions  become  singular  in  finite  time.  In  other  words,  the  shock  wave  solutions 
are  not  eigenfunctions,  they  are  states  of  degeneracy  of  the  eigenfunctions. 


However  with  the  introduction  of  the  forcing  function  periodic,  albeit  discontinu¬ 
ous,  solutions  may  be  found.  With  non-zero  forcing  function  (2.37)  has  the  following 
solution  for  the  first  mode, 


/(x)  =  'l7±^vfv/62  +  cos24 


where 


62=^{ 

32  l 


2  4u?  16 

r*  -I - L - 

9tt2  3tt2 


However,  for  mass  conservation  f(x)  must  satisfy 


f2w 

/  f(x)dx  =  0 
Jo 


(2.38) 


(2.39) 


(2.40a) 


w1±7V^(1  +  62)/0 


1  -  k2  sin2  xdx  =  0 


(2.406) 


where  k2  =  (1  +  62)  *.  (2.40b)  contains  an  elliptic  integral  of  the  second  kind  with 
modulus  k  which  has  a  real  solution  for  0  <  k2  <  1  only.  Therefore  there  are  no 


continuous  solutions  for 


12  ^2 

W<7V3 


(2.41) 


But  this  is  in  fact  the  region  of  interest;  the  neighborhood  of  resonance.  However,  as 
shown  by  Chester,  periodic  discontinuous  solutions  can  be  found  which  satisfy  (2.40). 
There  are  an  infinity  of  discontinuous  solutions  which  may  be  constructed,  using  the 
two  solutions  (2.38),  which  satisfy  (2.40).  However,  as  pointed  out  by  Chester,  when 
62  /  0  such  solutions  would  include  shocks  of  rarefaction  type,  which  are  non  'hysical. 
One  can  also  use  an  energy  argument  requiring  u  and  rj  to  increase  across  the  shock. 
Elaboration  on  these  arguments  are  given  in  Chester  [2],  and  Collins  5;,  in  his  more  gen¬ 
eral  study  of  nonlinear  wave  equations,  constructs  other  arguments  for  finding  unique 
discontinuous  solutions  based  on  an  “entropy”  argument  as  well  as  looking  for  discon¬ 
tinuous  solutions  as  limits  of  solutions  with  small  viscosity. 


Therefore,  taking  b2  =  0,  the  result  is  the  following  family  of  discontinuous  solu¬ 


tions, 


where 


/lx)  =  iVl(-^+cos  4] 


(2.42) 


sin"1  f-^1  <  *  <T  +  sin-1[-^l  (2.43) 

4\/6)  4  4y/6i 

for  -1  <  -  <  1.  Substituting  this  expression  into  (2.32)  results  in  the  first  order 

wave  height, 

4u\  4  [2 .  £ 

"'  =  ^-^nicosi  +  cos  4)  ,244> 

In  Figure  2.8  this  expression  is  plotted  for  fixed  x  (the  left  tank  wall)  as  a  function  of 
time.  The  forcing  function  is  equal  to  the  first  natural  frequency  (wi  =  0).  Since  the 
depth  has  been  ruled  out  through  the  limit  process,  all  one  can  say  is  that  the  depth 
is  “very  shallow”.  Chester  [2]  has  plotted  the  result  for  other  forcing  frequencies  by 
varying  wj.  They  are  similar  to  Figure  2.8  but  the  solutions  become  continuous  as 
exceeds  the  resonance  region,  given  by  (2.41). 

In  the  next  section  these  theories  will  be  compared  with  an  experimental  study. 
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3.  Experimental  Study  of  Elliptic  Standing  Waves  and  Hydraulic  Jumps 

3.1  Introduction 

An  experimental  study  was  undertaken  with  the  primary  goal  of  exploring  the 
nature  of  the  surface  waves  between  the  two  regions  discussed  in  the  previous  section. 
In  addition  some  previously  unrecognized  facets  of  elliptic  standing  waves  and  travelling 
hydraulic  were  illuminated.  The  experimental  observations  made  for  the  neighboorhood 
of  the  6  ~  t  plane  between  elliptic  standing  waves  and  shocks  will  be  discussed  in  detail 
in  Section  5.  It  is  the  purpose  here  to  describe  some  experimental  observations  obtained 
for  the  theoretical  regions  discussed  in  Section  2. 

The  regimes  observed  in  the  experiment  are  broken  down  in  Figure  3.1  where  a 
rough  diagram  of  the  waves  observed  for  various  combinations  of  e  and  6  are  plotted. 
The  points  near  the  diagonal  are  the  observed  values  of  t,  for  fixed  6 ,  where  standing 
waves,  of  any  variety,  were  observed  to  break  down.  Curves  of  e  =  6,  t  =  S2  as 
well  as  t  =  (n6)2  are  superimposed  on  the  diagram.  It  is  seen  that  the  breakdown  of 
standing  waves  occurs  at  e  ~  6  or  equally  well  c  —  (7r6)2.  There  is  then  an  intermediate 
region  where  compound  waves  are  observed,  which  will  be  discussed  in  Section  5,  and 
then  when  t  >  about  |  and  6  <  about  .15  the  travelling  hydraulic  jump  is  observed 
exclusively.  Although  the  travelling  jump  is  three  dimensionally  unsteady  in  a  strict 
sense,  in  the  large  it  is  essentially  two  dimensional  and  in  fact  very  stable.  Once  the 
hydraulic  jump  is  formed,  regardless  of  how  much  the  amplitude  of  the  tank  motion  is 
increased,  the  jump  is  maintained  with  the  strength  increasing  slightly  with  increasing 
tank  motion.  The  curved  labelled  line  in  the  figure  is  an  extrapolation  of  the  result  of 
Schwartz  and  Whitney  [14].  They  showed  that  for  fluid  of  infinite  depth  (the  purest  of 
elliptic  standing  waves)  the  maximum  height  of  a  stable  waveform  was  reached  when 
H 12 a  ~  .2.  In  terms  of  t  and  6  this  is  t  ~  .2/6  which  is  superimposed  on  the  figure. 
This  gives  an  upper  bound  for  elliptic  standing  waves,  and  may  be  an  upper  bound  for 
all  types  of  irrotational  waves  appearing  in  rectangular  vessels.  The  intersection  of  this 


curve  with  the  t  —  6 2  curve  at  6  ~  .6  suggests  that  for  6  >  .6  the  only  class  of  existing 
waves  are  the  elliptic  standing  waves. 


3.2  Experimental  Setup 


The  experimental  setup  is  depicted  in  Figures  3.2  and  3.3.  It  was  constructed  and 
setup  at  the  MIPAC  Research  Facility  at  the  Mathematics  Research  Center.  The  tank, 
constructed  of  thick  plexiglass,  has  the  following  dimensions:  tank  length,  2a  =  18", 
tank  breadth  =  3",  and  the  tank  height  =  15".  Although  the  tank  height  was  18",  the 
fill  depth  was  limited  to  about  6"  due  to  the  large  inertial  force  imposed  by  the  fluid  on 
the  driving  mechanism  causing  perturbations  to  the  otherwise  strictly  sinusoidal  forcing 
function.  The  natural  frequencies  in  Hz  of  the  first  five  (two  dimensional)  modes,  as 
a  function  of  the  dimensionless  depth  6,  for  a  vessel  of  this  dimension  is  shown  in 
Figure  3.4.  It  is  important  to  note,  however,  that  in  a  rigid  tank  undergoing  forced 
harmonic  oscillation  only  odd  modes  may  be  excited.  In  fact,  it  appears  that  the  only 
way  to  excite  even  modes  is  to  move  the  endwalls  harmonically.  This  is  how  Taylor 
[  16l  excited  standing  waves  in  his  study  of  the  limiting  configuration  of  standing  waves 
in  deep  water.  In  his  case  the  mode  studied  and  depicted  in  photos  in  his  paper  was 
the  m=2  mode.  Due  to  the  rigid  walls  only  odd  modes  are  elicited  in  this  study. 

A  steel  shaft,  rigidly  bolted  to  the  bottom  center  of  the  tank  and  mounted  in 
rotary  bearings  on  a  wooden  support  stand,  acts  as  the  pivot  for  the  roll  motion  of 
the  tank.  A  driving  arm,  attached  perpendicular  to  the  shaft,  is  connected  to  the 
driving  mechanism,  which  oscillates  vertically.  With  this  arrangement  the  vertical 
oscillation  of  the  driving  mechanism  imparts  an  angular  rotation  to  the  tank,  i.e.  roll 
motion.  The  driving  mechanism  is  a  Bruel  &  Kjaer  Vibration  Exciter  Type  4808.  It  is 
an  electromagnetic  vibration  exciter  much  like  a  loudspeaker,  where  the  movement  is 
produced  by  a  current  passing  through  a  coil  in  a  magnetic  field.  The  amplitude  of  the 
excitation  is  varied  by  controlling  the  current.  The  harmonic  output,  to  the  driving 
mechanism,  is  created  in  and  controlled  by  the  HP  5451C  Fourier  Analyzer  shown  in 
Figure  3.3. 
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Figure  3.2  The  apparatus  used  to  study 
the  waves  in  a  rectangular  basin.  The 
black  object  in  the  background  is  a  ver¬ 
tical  oscillator  which  provides  the  forc¬ 
ing  function.  The  electronics  to  the  left 
pre-process  the  wave  height  time  series 
for  input  into  the  Fourier  analyzer. 


Figure  3.3  The  HP  Fourier  analyzer 
which  drives  the  oscillator;  records 
and  analyzes  the  experimental  data. 


The  HP  5451C  Fourier  Analyzer  is  a  programmable  minicomputer  which  creates, 
inputs,  outputs,  and  analyzes  time  series.  Fourier  analysis  programs  such  as  the  (a) 
forward  and  inverse  Fourier  transform,  (b)  power  and  cross-power  spectrum,  (c)  trans¬ 
fer  function,  (d)  coherence  function,  (e)  convolution,  (f)  auto  and  cross  correlation,  (h) 
ensemble  averaging,  etc.,  are  built  in  and  can  be  performed  through  keyboard  com¬ 
mands  without  special  software.  In  addition,  the  analyzer  can  be  programmed  to  run 
the  experiment  and  input  and  analyze  the  various  sensory  measurements  such  as  the 
local  wave  height  and  the  tank  rotary  displacement.  In  this  experiment  the  Fourier 
analyzer  is  programmed  to  do  the  following.  It  creates,  given  the  desired  frequency 
or  frequencies,  and  outputs  a  time  series  to  the  exciter.  The  amplitude  of  the  exciter 
is  controlled  in  real  time  and  can  be  continuously  varied.  The  Fourier  analyzer  then 
recieves  and  analyzes  on  command  the  time  series  of  the  wave  height  and  or  tank  dis¬ 
placement.  A  hard  wired  external  plotter  is  used  to  plot  the  results,  and  an  hard  disk 
is  used  for  storage  of  the  programs  and  data  blocks. 

The  wave  height  is  measured  using  a  taunt-wire  capacitive  wave  guage  mounted  one 
inch  from  the  left  tank  wall.  The  electronics,  which  convert  the  changing  capacitance, 
produced  by  the  water  riding  up  and  down  the  wire,  to  a  changing  voltage,  which  is 
then  input  into  the  Fourier  analyzer,  is  mounted  to  the  left  of  the  tank  in  Figure  3.2. 
The  Fourier  analyzer  plots  the  signal  and  computes  the  Fourier  transform  to  deduce 
the  frequency  content. 


3.3  Experimental  Observations 

In  the  experiment  the  wave  forms  are  elicited  by  forcing  the  tank  to  oscillate. 
Therefore  there  is  not  complete  congruence  between  the  solutions  observed  and  the 
bifurcation  analysis  of,  for  example,  Tadjbakhsh  and  Keller.  However,  when  the  forc¬ 
ing  frequency  is  near  resonance,  the  solutions  of  the  forced  problem  should  be  very 
near  those  of  the  bifurcation  analysis.  Therefore  the  forcing  frequency  will  be  near 
resonance  unless  otherwise  noted.  An  experimental  setup  of  this  kind  is  dependent  on 
five  dimensionless  parameters;  7  and  d  which  specify  the  point  of  rotation  (however  in 
this  experiment  these  are  fixed  at  7  =  |  and  d  =  6),  6  the  fluid  depth,  and  u  and  t? 
the  frequency  and  amplitude  of  the  forcing  function.  The  amplitude  of  the  response, 
c,  is  also  a  parameter,  but  it  is  posited  that  there  is  a  functional  relationship  between 
e  and  t?.  Therefore  for  a  particular  “run”,  once  the  depth  is  fixed  the  only  parameters 
are  the  frequency  u>  (actually  /  will  be  used  where  /  =  u;/ 2n)  and  amplitude  t?  of  the 
forcing  function. 

Observations  of  elliptic  standing  waves  will  now  be  discussed.  As  pointed  out  in 
Section  2  these  solutions  occur  at  all  depths  but  for  only  a  limited  range  of  amplitude. 
Therefore  consider  the  fill  depth  6  =  .2.  Here  the  relevant  natural  frequencies  are 
fi  ~  1.04 Hz,  /3  =  2.22 Hz,  and  /s  =  2.93 Hz.  In  Figures  3.5,  3.6,  and  3.7  photos  of 
the  1,  3,  and  5  mode  respectively  are  shown.  In  each  case  the  forcing  function  was 
set  equal  to  the  linear  natural  frequency.  These  figures  may  be  compared  to  the  result 
of  the  theory  of  Tadjbakhsh  Sc  Keller  shown  in  Figures  2.4  and  2.6.  The  results  are 
similar.  These  figures  are  in  fact  typical  of  the  elliptic  variety  of  standing  waves.  It  is 
interesting  that  the  waves  of  this  type  look  elliptic.  The  body  of  fluid  moves  uniformily 
as  if  every  particle  were  aware  of  every  other  particles  motion.  A  striking  contrast  to 
hyperbolic  waves. 


Figure  3.5  A  classical  mode  1  elliptic  standing  wave,  where  b  -  .20,  /  -  / i 


Figure  3.6  A  mode  3  elliptic  standing  wave,  showing  signs  of  a  three  dimen 
sional  instability  at  the  peaks,  with  6  -  .2  and  /  f3  2.22  Hz. 


Figure  3.9a  With  6  -  .077  and  /  =  /5  =  2.68  a  mode  5  elliptic  standing  wave 
occurs. 


Figure  3.9h  However,  at  the  same  depth  6  .077.  but  with  a  lower  frequency 

/  /,  .64  and  sufficiently  large  amplitude,  a  travelling  hydraulic  jump  is 

formed. 


The  time  series  for  the  wave  height  shows  this  as  well.  Figure  3.11a  shows  the  time 
series  for  the  “deep  water  jump”  when  6  =  .104  and  f  —  J\.  Although  the  time  series 
is  “almost  discontinuous”,  and  periodic  “in  the  large”  there  is  some  irregularity  due  to 
the  small  amount  of  dispersion  present.  The  Log  of  the  Fourier  transform  of  this  wave 
height  is  given  in  Figure  3.1  lb.  It  shows  the  richness  of  the  harmonics  present.  The  time 
series  for  the  “shallow  water  jump”  is  depicted  in  Figure  3.12a.  Here  the  time  series  is 
more  regular  and  the  “almost  discontinuous”  nature  of  the  wave  height  is  obvious.  The 
presence  of  dissipation  prevents  the  solution  from  becoming  truly  discontinuous  but 
nevertheless  the  experimental  time  series  is  approaching  the  theoretical  prediction  of 
Chesters  shown  in  Figure  2.8.  The  Log  of  the  Fourier  transform  for  this  time  series  is 
shown  in  Figure  3.12b.  The  explosion  of  harmonics  over  that  occurring  for  the  standing 
wave  problem  is  obvious.  But  this  is  to  be  expected  from  a  function  which  is  almost 
discontinuous. 

In  the  next  section  a  theory  for  the  “missing  link”  between  elliptic  standing  waves 
and  the  hydraulic  jump  is  proposed. 
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Figure  3.11a  Time  series  for  a  “deep  water”  travelling  hydraulic  jump  where 
b  -  .105  and  /  =  /,  =  .74.  At  this  depth  there  remains  enough  dispersion 
to  produce  an  irregular  travelling  jump. 
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Figure  3.11b  Log  magnitude  of  the  Fourier  transform  of  the  time  series  in 
Figure  3.11a.  Due  to  the  “almost  discontinuous”  nature  of  the  time  series 
the  spectrum  is  rich. 
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Figure  3.12a  Time  series  for  a  “shallow  water”  travelling  hydraulic  jump 
where  6  =  .077  and  f  -  f\  -  .64.  Here  the  function  is  more  regular  and 
begins  to  resemble  the  shallow  water  solution  of  Chester  in  Figure  2.8. 
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figure  3.12b  Log  magnitude  of  the  Fourier  transform  of  the  time  series  in 
Figure  3.11b.  The  spectrum  is  again  rich  due  to  the  nature  of  the  hydraulic 


4.  Theory  -  II:  Cnoidal  Standing  Waves 


The  solution  for  finite  amplitude  standing  waves  in  fluid  of  arbitrary  but  finite 
depth  was  first  derived  by  Tajdbakhsh  and  Keller  [15]  and  extended  to  arbitrary  mode 
number  in  Section  2.  It  was  also  shown  in  Section  2  that  this  solution  is  not  uniformily 
valid  in  the  amplitude.  When  e  =  0(62)  there  is  a  spillover  of  the  higher  harmonics 
onto  the  zeroth  order  solution  violating  the  tenet  of  the  expansion.  This  spillover  is 
in  fact  the  first  sign  of  the  change  of  type  of  the  standing  wave.  The  amplitude  (e) 
is  beginning  to  get  large  enough  to  challenge  the  role  of  the  dispersion  (62).  When 
t  <<  6  2  the  dispersion  dominates  and  the  solution  is  the  well-behaved  elliptic  standing 
wave,  and  when  t  >>  62  the  dispersion  can  no  longer  check  the  coalescence  of  the 
harmonics  and  they  run  together  as  a  hydraulic  jump. 

There  is  however  an  intermediate  region  where  the  dispersion  and  the  amplitude 
are  in  harmony,  i.e.  e  =  0'62).  A  theory  for  the  wave  forms  in  this  region  is  proposed 
in  this  section.  In  order  for  the  relation  f  =  0(6 2)  to  hold  6  must  be  sufficiently  small. 
Therefore  an  asymptotic  solution  is  sought  with  the  simultaneous  limit  e  —*  0  and 
6  — *  0,  but  with  the  requirement  that  e/62  remain  constant.  This  is  done  by  taking 
the  relationship  6  -  v/re  where  r  =  0(1)  and  carries  the  sign  of  e. 

The  governing  equations  and  boundary  conditions  are  the  same  as  (2.l)-(2.3)  but 
with  the  substitution  6  =  \Zre, 
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With  f  as  the  only  small  parameter  appearing,  a  power  series  solution  in  e  is 
sought, 
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V  =  Vo  +  (-V1  +  <2f?2  +  ••• 
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(4.5) 

(4.6) 
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These  expressions  are  substituted  into  the  governing  equations  and  boundary  condi¬ 
tions,  the  free  surface  boundary  conditions  are  expanded  in  a  Taylor  series  about  the 
still  water  level,  and  terms  proportional  to  like  powers  of  t  are  set  to  zero  resulting  in 
a  sequence  of  boundary  value  problems.  The  Laplace  equation  gives, 
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etc..  These  may  then  be  substituted  into  the  free  surface  boundary  conditions.  The 
zeroth  order  problem  is  then 
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and  this  may  be  continued  to  higher  order.  Here  the  solution  will  be  carried  through 
the  first  order. 


Combining  (4.11)  and  (4.12)  the  zeroth  order  problem  is  a  homogeneous  wave 
equation, 


U0o  =  0 


(4.15) 


where  U  is  the  D’Alembertian, 

2  &  d2 

U~UJ°dt 2  dx2 


(4.16) 


with  the  additional  requirement  that  d^o/dx  vanish  on  the  vertical  boundary  x  —  ±  \. 
(4.15)  and  the  wall  boundary  condition  have  the  general  solution 


uq  =  am  (4.17) 

0o  =  fit)  +  fix)  (418) 


where  £  =  t  +  amx,  \  -  t  -  amx,  am  =  mir,  and  x  =  x  +  With  (4.11)  the  zeroth 
order  wave  height  is  found  to  be 


rio  =  -am[/'(0  +  /'(*)] 


(4.19) 


The  function  /(£)  is,  however,  unknown  at  this  order.  The  Fredholm  solvability  con¬ 
dition  applied  at  the  next  order  will  result  in  an  expression  for  /(£). 

The  zeroth  order  result  (4.17)-(4.18)  are  substituted  into  (4.13)-(4.14).  The  two 
equations  are  combined  to  form  the  following  inhomogeneous  wave  equation  for  0 1, 


U0i  =  F\(x,t) 


(4.20) 


where 

FliX't)  =  ^j=  {/""({)  +  !""{x)}  -  2amWl{/"(0  +  /"(X)} 
-  3a»,{/'(e)He)  +  /'(x)/"(x)} 

+  «3n,{/"(o/'(x)  +  no/w(x)} 


(4.21) 


(4.20)  has  a  solution  if  and  only  if  Fi(x,t)  satisfies, 
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d£  =  0 


(4.22) 


This  is  a  variant  of  the  Fredholm  alternative  and  is  proved  in  Appendix  III.  Application 
of  this  requirement  results  in  the  following  nonlinear  ordinary  differential  equation  for 


^/""(X)  -  2amWl/'#(x)  -  2a3m/'(X)/"(x)  =  0 

but  defining  h(\)  =  -am/'(x)  and  integrating  the  equation  once, 

rh"(x)  "  £«x)  +  =^-*2(x)  -  jK-h  =  0 

am  2am  4am 


(4.23) 


(4.24) 


Conservation  of  mass  requires  h(x)dx  —  0,  therefore  I\  is  a  positive  definite 

constant  given  by 

r  20tm 

I\  —  I  h2(x)dx  (4-25) 

Jo 

It  is  important  to  note  that  the  nonlinearity  of  the  resulting  expression  for  h(\)  is 
quadratic.  This  suggests  that  the  nonlinearity  of  the  operator  is  quadratic  in  that 
region  of  the  6  -  e  plane  where  €  =  0(6 2 )  (and  this  may  also  be  true  for  c  >>  62).  The 
operator  for  f  <<  62  is  cubic  in  the  nonlinearity  as  discussed  in  Section  2. 

The  solution  of  (4.24)  can  be  found  in  terms  of  Jacobian  elliptic  functions,  i.e., 
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h{x)  =  A  +  £cn2(x;*) 


(4.26) 


when  u’i  =  |omA  -  |ra^,(l  -  2«2),  and 


A  =  ~\a2mr{K2  -  1  +  -E(tt;k)} 
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and  the  constant  1\  is  found  to  have  the  following  expression, 
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The  cn(x;«)  is  the  family  of  cnoidal  functions  with  modulus  k.  Necessary  details  for 
this  and  the  other  Jacobian  elliptic  functions  are  given  in  Appendix  IV.  f{\)  is  found 
by  integrating  h(x) 

4 

/(x)  =  — TamK2Z(x'iK)  +  constant  (4.30) 

3 

where  Z(x; «)  is  the  Jacobian  Zeta  function.  The  constant  is  arbitrary  and  it  is  chosen 
for  convenience  to  be  zero  to  satisfy  /02“m  f{x)dx  =  0.  A  complete  description  of  the 
zeroth  order  solution  may  now  be  given, 


w0  =  (4.31) 

4 

4>o(x,  t)  =  -  -t amK2\Z(t  +  amx\  k)  +  Z(t  -  amx;  k)]  (4.32) 
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9o(*,0  =  --o^r[K2  -  1  +  -E(x;ie)] 
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+  ^K2raj^[cn2(t  +  amx;  k)  +  cn 2(t  -  amx; «)]  (4.33) 

O 

and  the  first  order  correction  to  the  frequency  is  found  to  be 

w,  =  -2ra^[^E(^;«)  -  J(2  - /c2)]  (4.34) 

7T  6 


which  is  positive  definite.  Therefore  the  bifurcation  is  supercritical;  i.e.  the  natural 
frequency  increases  with  increasing  amplitude.  It  is  interesting  to  note  that  t  and 
usually  appear  together.  Recalling  that  c  =  |62,  it  was  found  in  the  solution  for 
elliptic  standing  waves  that  spillover  occurs  when  e  =  (om6)2.  In  addition,  as  shown  in 
Appendix  IV  equation  (IV.6),  the  cnoidal  function  has  a  complete  Fourier  cosine  series 
showing  that  when  i  =  0(62)  all  the  higher  harmonics  are  present  in  the  solution. 

With  /(x)  known,  it  may  be  substituted  into  the  right  hand  side  of  (4.20)  and 
upon  integration  the  first  order  solution  for  is  found  to  be 

=  ■  X)  +  ^-imnx)  +  Hx)/'{() ]  +  |/.({)  +  /,<x)l 


(4.35) 


and  by  combining  with  (4.9)  the  complete  expression  for  the  first  order  potential  is 


*.  =  -  ~rU  +  x)  +  ^1  /(£)/'<*)  +  /(x)f'U)  I  +  +  s/)2[fc'(e)  +  fc'(x)l 

+l/i(£)  +  /i(x)l  (4-36) 

and  with  (4.13)  the  first  order  wave  height  is 

1,  =  +  — |M£)  +  fe(x)]  -  ’M£)Mx)  -  =»V(0  +  *"(*)] 

2am  am  22 

+  ~[/(«fc'(x)  +  /(x)fc'(£)]  -  ||MC)  -  Mx)P 

-«m[/!(£)  +  /!(x)l  (4.37) 

The  function  / j(£),  as  well  as  the  second  order  correction  to  the  frequency,  may  be 
found  by  imposing  the  Fredholm  solvability  condition  to  the  second  order  problem. 

The  main  result  is  that  the  zeroth  order  wave  height  is  given  by  a  family  of  left 
and  right  running  cnoidal  waves  with  parameter  k  varying  over  0  <  k2  <  1.  The 
combination  of  the  left  running  and  right  running  waves  forms  a  standing  cnoidal  wave. 
If  k 2  =  1,  then  cn(£;  1)  — »  sech(£)  which  would  then  suggest  that  the  waves  would  be 
a  combination  of  a  right  and  left  running  solitary  wave.  However,  this  is  not  the  case 
since  the  solitary  wave,  due  to  the  necessity  of  an  infinitely  long  tail,  will  not  satisfy 
the  periodicity  condition.  In  fact  if  it  is  required  that  the  solution  be  periodic  in  time 
with  period  2t\  then  there  is  only  one  member  of  the  k  family  which  will  satisfy  this. 
The  cn2(£; k)  function  is  periodic  with  period  2K( |; «),  where  K(|;«c)  is  the  complete 
elliptic  integral  of  the  first  kind.  Therefore  the  admissable  value  of  K(|;  k)  is  given  by 
the  root  of  the  transcendental  equation, 

/2  =0  (4-38) 

Jo  v  1  -  sin  c 

A  numerical  solution  of  this  equation  results  in  a  value  for  k 2  of  k 2  ~  .9691. 
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The  second  order  frequency  correction  u>2  and  the  function  /i(x)  may  be  found 
through  application  of  the  solvability  condition  on  the  second  order  problem.  The 
second  order  problem  has  the  form 

d<t> i  ,  d<j> o  ,  dfodfi  1  / d<f>i\ 2  d2<f>i  ,  x 

^  +  ..0-^+^—  +—  —  +  Tr(-g~)  +^0^=0  (4.39) 


dr)2  ld(f>z  dr)  i  dr\0  dr)0  d<t>\  dt\ \d<j>o 

<jJ°  dt  t  dy  Wl  dt  + W2  dt  +  dx  dx  dx  dx 

_  rn  _  1i  d2<t>  i  _  0 

T  dy 2  T  dy2 


(4.40) 


When  these  equations  are  combined  and  the  solvability  condition  is  applied,  after  some 
integration  and  algebraic  manipulation  the  resulting  differential  equation  for  /i(x)  >s 

r/r(x)  -  ^/((x)  +  ^2“Mx)/((x)  =  constant 
a*  2a£, 


.6W2  3/j  4rdi 

+  / - ±  -l.  - -  -i - 

am  I5am 


}MX)  -  -  ~h'(X)h'(X) 


(4.41) 


where  the  numbers  dx  and  d2  are  given  in  Appendix  IV,  and  the  integration  constant 
is  chosen  such  that  the  integral  of  /i(x)  is  zero.  The  equation  is  a  linear  differential 
equation  with  nonconstant  coefficients  which  has  the  homogeneous  solution  /{/,(x)  — 
fih(x)  where  /?  is  an  arbitrary  constant.  The  inhomogeneous  part  must  satisfy  an 
additional  solvability  condition.  After  multiplication  of  (4.41)  by  h(x)  and  integration 
from  0  to  2am  the  left  hand  side  is  zero  and  the  integral  of  the  right  hand  side  results 
in  an  expression  for  u>2, 
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+  (5^2/4  -  2/5) 


(4.42) 


where  the  constants  d  1,  d2,  I4,  and  Is  are  given  in  Appendix  IV.  After  substitution  of 
this  into  (4.41)  it  can  be  shown  that  (4.41)  has  the  general  solution 
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The  arbitrariness  of  Aj  is  due  to  the  integration  constant  in  (4.41)  and  the  arbitrariness 
of  .42  is  due  to  the  homogeneous  solution  of  (4.41).  This  arbitrariness,  however,  may 
be  eliminated.  To  satisfy  conservation  of  mass  the  integral  of  f[{x)  in  0  to  2am  must 
be  zero,  therefore 
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(4.44) 


and  if  c  is  defined  such  that  it  is  proprotional  to  the  leading  order  wave  height  then 
the  inner  product  of  »/o  with  f/j  should  be  zero.  This  gives  a  value  for  A?,  of 


‘ham 


(4.45) 


Therefore  substitution  of  (4.43)-(4.45)  into  (4.36)  and  (4.37)  gives  a  complete  descrip¬ 
tion  for  the  first  order  solution. 

In  Figure  4.1  a  bifurcation  diagram  for  the  frequency  is  shown,  cjj  is  positive 
definite  and  it  appears  that  cj2  is  as  well  so  the  bifurcation  is  subcritical  for  all  the 
modes;  the  natural  frequency  of  the  cnoidal  standing  waves  increases  with  increasing 
amplitude.  A  typical  time  series  for  the  wave  height  at  the  left  tank  wall  is  given  in 
Figure  4.2. 


In  Figure  4.3,  the  evolution  of  the  mode=l  cnoidal  standing  wave  is  shown.  The 
result  is  plotted  as  a  function  of  space  with  each  successive  frame  of  a  temporal 
period.  This  wave  is  very  similar  to  that  found  by  Vanden-broeck  and  Schwartz  [18], 
using  a  numerical  technique,  and  shown  in  their  Figure  3.  Figures  4.4  and  4.5  show 
the  evolution  of  a  2-mode  and  3-mode  respectively.  The  “travelling  character”  of  these 
waves  is  to  be  contrasted  with  the  “standing  character”  of  the  elliptic  standing  waves. 
In  Figure  4.6  the  evolution  of  a  1  mode  elliptic  standing  wave  is  plotted.  Here  the 
wave  appears  to  move  vertically  and  otherwise  remain  coherent.  Figure  4.7  shows  the 
evolution  of  a  3-mode  elliptic  standing  wave. 


Figure  4.1  Bifurcation  diagram  for  the  natural  frequency  (the  first  five  modes) 
of  standing  cnoidal  waves. 


Figure  4.2  Time  series  for  the  wave  height  at  the  left  tank  wall  for  a  cnoidal 
standing  wave. 


t  =  o 


Figure  4.6  Temporal  evolution  of  the  x- 
distribution  of  the  wave  height  for  a  1- 
mode  elliptic  standing  wave. 


t  =  7T 


The  distribution  of  the  2-mode  at  t  =  ir  at  the  bottom  of  Figure  4.4  raises  an 
interesting  question  about  the  enclosed  angle  of  the  crest  of  a  large  amplitude  cnoidal 
standing  wave.  It  is  well  known  that  the  Stokes  wave  has  a  limiting  angle  of  120  . 
Schwartz  and  Whitney  [14]  have  shown  that  the  elliptic  standing  wave  probably  has  a 
limiting  angle  of  90°.  The  cnoidal  standing  wave  appears  to  have  a  limiting  angle  less 
than  90°,  probably  much  less  than  90°,  and  possibly  a  limiting  configuration  that  is  a 
cusp. 

The  cnoidal  standing  wave  is  the  periodic  surface  wave  occurring  when  the  fluid 
is  shallow  and  £  =  0(62).  When  c  >>  62  it  was  shown  that  no  periodic  solutions  exist 
to  the  innate  problem  although  a  hydraulic  jump  is  formed  in  response  to  a  forcing 
function.  There  is  obviously  a  vague  region  remaining  between  e  —  0(62)  and  t  >>  62. 
The  transition  from  a  cnoidal  standing  wave  to  a  hydraulic  jump  is  an  interesting 
but  still  unclear  theoretical  problem.  In  the  next  section  experimental  observations  of 
cnoidal  standing  waves  as  well  as  another  class  of  waves  witnessed  between  £  =  0(S2) 
and  £  >>  S2  will  be  elaborated. 


5.  Experimental  Study  of  the  Transition  Region 

The  transition  region  is  defined  to  be  the  region  in  the  amplitude,  t ,  with  fixed  and 
small  6,  after  the  breakdown  of  the  elliptic  standing  waves  and  before  the  formation 
of  the  hydraulic  jump.  It  was  shown  in  the  previous  section  that  the  eigenfunction  for 
the  nonlinear  problem  in  this  region  is  the  cnoidal  standing  wave.  The  first  step  is  to 
see  if  such  waves  are  observed  in  the  experiment.  It  will  be  shown  that  this,  in  fact, 
is  the  case.  The  parameters  are  chosen  such  that  an  isolated  1-mode  cnoidal  standing 
wave  occurs.  The  depth  paramter  is  6  ~  0.2,  and  the  forcing  frequency  is  /  =  1.0 
which  is  just  above  the  first  linear  natural  frequency  f\  ~  0.97.  The  photos  in  Figures 
5.1,  5.2,  and  5.3  are  of  this  wave  at  three  different  times.  The  wave  has  the  distinct 
undulating  motion  associated  with  the  cnoidal  standing  wave  and  not  observed  in  the 
elliptic  standing  wave  which  has  (almost)  a  node  at  the  center  of  the  tank.  Comparing 
Figures  5. 1-5.3  with  Figure  4.3  it  is  seen  that  5.2  resembles  the  evolution  in  4.3  at 
t  -  (9  or  10)-^,  5.2  at  t  ~  (6  or  7)y^  and  5.3  at  <  =  f  •  Figure  5.4a  is  the  time  series 
for  the  wave  height  at  the  left  tank  wall  for  the  wave  in  Figure  5. 1-5.3  and  Figure  5.4b 
is  the  Log  of  the  Fourier  transform  of  this  time  series.  The  time  series  closely  resembles 
that  predicted  by  the  cnoidal  wave  theory  in  Figure  4.2. 

The  eigenfunction  analysis  of  the  previous  sections  has  suggested  that  the  eigen¬ 
functions  when  t  <<  62  are  the  elliptic  standing  waves,  and  when  t  =  0(62)  the 
eigenfunctions  are  the  cnoidal  standing  waves,  and  when  c  >>  b2  it  was  suggested  that 
there  are  no  periodic  characteristic  solutions.  It  is  possible  that  there  is  another  class 
of  eigenfunctions  between  e  =  0(62)  and  t  >>  62  but  this  is  unlikely.  These  remarks 
serve  as  an  introduction  to  the  experimental  observations  which  follow. 
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Figure  5.4a  The  time  series  for  the  wave  height  at  the  left  tank  wall  cor¬ 
responding  to  Figures  5. 1-5.3.  The  parameters  are  6  ~  0.2,  f\  =  .97,  and 
/  =  1.0. 


Figure  5.4b  The  Log  magnitude  of  the  Fourier  transform  of  the  time  series 
in  Figure  5.4a. 
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The  complex  waves  observed,  as  the  amplitude  increases  and  the  “clean”  1-mode 
cnoidal  standing  wave  changes  into  a  wave  field  with  greater  spatial  complexity,  are 
probably  due  to  the  selecting  process  of  the  forcing  function  which  is  present.  So  rather 
than  a  new  class  of  eigenfunctions  these  waves  are  probably  compound  cnoidal  waves 
and  in  fact  the  observations  lend  credance  to  this  judgement.  There  are  two  distinct 
compound  waves  which  occur  so  often  (for  a  large  but  compact  range  of  parameters) 
that  they  may  be  given  titles.  They  are  the  compound  1-3  wave  and  the  compound 
wiggle  wave.  The  compound  1-3  wave  appears  at  shallow  depths  but  not  too  shallow 
(6  >  0.15)  and  looks  like  a  cnoidal  1-mode  joined  with  a  cnoidal  3-mode.  Figures  5.5 
and  5.6  are  photos  of  the  compound  1-3  wave  at  two  different  times  showing  (almost) 
the  symmetry.  The  forcing  is  at  the  first  natural  frequency  but  the  presence  of  the 
third  mode  is  quite  obvious.  This  wave  has  a  rather  interesting  time  series  as  well  and 
it  is  shown  in  Figure  5.7a  with  its  Fourier  transform  in  Figure  5.7b.  The  time  series  is 
periodic  with  an  obvious  presence  of  higher  harmonics. 

The  compound  wiggle  wave  occurs  at  shallower  depths  at  sufficient  amplitude. 
Figures  5.8-5.11  show  an  example  of  this  wiggle  wave  at  four  different  times.  The 
depth  is  6  =  .084  and  the  forcing  frequency  equals  the  linear  natural  frequency.  The 
wave  is  periodic,  repeatable,  and  deterministic  but  has  many  spatial  harmonics  (most 
notably  a  5-mode  in  Figure  5.10)  and  the  harmonics  travel  at  different  speeds  from 
the  fundamental.  A  time  series  for  this  wiggle  wave  is  given  in  Figure  5.12a,  with  its 
Fourier  transform  in  Figure  5.12b,  where  the  “wiggles”  are  quite  apparent.  The  Fourier 
transform  also  suggests  that  these  harmonics  are  integral  multiples  of  the  fundamental. 
It  is  interesting  to  note  that  the  compound  wiggle  wave  is  always  a  prelude  to  the 
hydraulic  jump.  If  a  compound  wiggle  wave  is  formed  and  the  amplitude,  of  the 
forcing  function,  is  increased  further  a  threshhold  is  reached  where  the  wave  changes 
(bifurcates?)  into  a  travelling  hydraulic  jump. 
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Figure  5.7a  Time  series  at  the  left  tank  wall  for  the  compound  1-3  standing 
wave.  The  parameters  are  6  ~  0.17,  f\  —  .91,  and  /  =  .95. 
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Figure  5.7b  The  Log  magnitude  of  the  Fourier  transform  of  the  time  series 
in  Figure  5.7a. 


Figure  5.9  A  compound  wiggle  wave  with  the  same  parameters  as  Figure  5. 
at  a  different  time  showing  how  the  wave  crests  travel  independently. 
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Figure  5.12a  Time  series  for  the  wave  height  at  the  left  tank  wall  for  the 
compound  wiggle  wave.  The  parameters  are  6  =  .098,  fx  =  .716,  and  /  =  .67. 
The  wave  form  is  living  up  to  its  name  but  remains  periodic. 
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Figure  5.12b  The  Log  magnitude  of  the  Fourier  transform  of  the  time  series 
in  Figure  5.12a. 


Figure  5.13  shows  an  example  of  another  variety  of  compound  wiggle  wave  in 
deeper  water.  Figures  5.14-5.17  are  four  other  interesting  time  series  which  were  ob¬ 
served  and  are  worthy  of  mention.  Figures  5.14,  5.15  are  of  two  different  compound 
wiggle  waves  showing  the  increase  in  harmonic  complexity  as  the  amplitude  increases 
and  Figure  5.16,5.17  are  time  series  of  the  wave  height  for  ordinary  1-mode  cnoidal 
standing  waves  as  they  become  more  nonlinear.  These  observations  lead  to  a  cascade 
hypothesis. 

The  following  scenario  may  be  constructed  based  upon  the  observations  recorded 
in  the  figures  in  this  section.  When  c  =  0(52)  the  solution  observed  is  that  of  a  cnoidal 
standing  wave  which  closely  resembles  the  eigenfunction.  As  the  amplitude  is  increased 
such  that  it  is  slightly  greater  than  0(5 2)  the  third  mode  natural  frequency  becomes 
closer  to  three  times  the  first  mode.  This  brings  in  the  3-mode  cnoidal  standing  wave 
forming  the  compound  1-3  standing  wave  observed  in  Figures  5.5-5. 7.  Now  as  the 
amplitude  is  increased  further  relative  to  5 2  a  point  is  reached  such  that  the  5-mode 
cnoidal  standing  wave  is  excited.  This  produces  the  compound  wiggle  wave  observed 
in  Figures  5.7-5.14.  As  the  amplitude  is  further  increased  relative  to  5 2  the  cascade  of 
higher  harmonics  continues  until  the  dispersion  is  so  small  that  it  no  longer  can  check 
the  flood  of  harmonics  and  a  travelling  hydraulic  jump  is  formed. 

In  other  words  as  c  increases  to  e  >>  52  with  52  sufficiently  small,  the  higher  natu¬ 
ral  frequencies  become  closer  and  closer  to  being  integral  multiples  of  the  fundamental, 
and  apparently  this  occurs  in  order  of  increasing  mode  number  causing  the  cascade 
process  described  above. 


Figure  5.14  A  variation  of  the  time  series  for  the  compound  wiggle  wave. 
The  parameters  are  6  =  .08,  /j  =  .6655,  and  /  =  .67.  This  time  series 
corresponds  to  the  photos  in  Figures  5.8-5.11. 


0  1  2  3  4  5 

Time,  sec 


Figure  5.15  The  time  series  for  a  further  variation  of  the  compound  wiggle 
wave  when  6  -  0.07,  /j  =  .61,  and  /  =  .6.  Although  the  wave  remains 
periodic  there  is  a  rich  set  of  harmonics  present. 
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Figure  5.16  A  time  series  for  a  1-mode  cnoidal  standing  wave  of  large  am¬ 
plitude.  Although  (spatially)  in  appearance  this  wave  resembles  a  1-mode 
cnoidal  standing  wave  it  has  reached  sufficient  amplitude  to  require  higher 
order  terms.  The  parameters  are  6  =  .14,  =  .84,  and  /  =  .84. 
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Figure  5.17  A  time  series  for  another  variation  of  the  nonlinear  1-mode 
cnoidal  standing  wave.  The  parameters  are  6  =  .112,  =  .762,  and  /  =  .8. 
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Appendix  I 

A  coordinate  system  S(X,Y,Z)  is  fixed  in  Newtonian  space  and  a  second  coor¬ 
dinate  system  S'(x,y,z)  is  attached  to  the  center  of  the  undisturbed  water  level  of  a 
three  dimensional  partially  filled  cuboidal  container.  From  the  theory  of  dynamics  it 
is  well  known  that  when  the  coordinate  system  S'  is  accelerating  relative  to  the  fixed 
system  S,  the  relationship  between  the  absolute  and  relative  acceleration  is 

A  =  a  +  2nxu  +  nxr  +  rixnxr  +  nxd  +  rixnxd  (/. 1) 

where  O  is  the  angular  velocity  of  the  vessel  and  d  is  the  fixed  distance  between  the 
origins  of  S  and  S'.  This  expression  is  substituted  into  the  Euler  equations  of  motion 
for  an  inviscid  fluid  to  obtain  the  governing  equations  for  the  fluid  in  the  vessel  relative 
to  S', 

^  +  2nxu  +  nx(r  +  d)  +  nxftx(r  +  d)  = 

~-Vp+gVh  (/.  2) 

P 

where  h  is  the  elevation  above  some  datum.  Because  of  the  presence  of  fl  the  fluid 
motion  is  not  irrotational.  However  if  it  is  assume^  that  the  vorticity  is  strictly  due  to 
the  rigid  body  motion  of  the  vessel  the  velocity  field  can  be  decomposed  into 

u  -  V<*>-  rixr  (7.3) 

and  the  motion  relative  to  S'  is  assumed  to  be  irrotational.  The  governing  equation  is 
then  the  Laplace  equation  in  the  interior  of  the  fluid  and  substitution  of  (1.3)  into  (1.2) 
results  in  a  modified  Bernoulli  equation  for  the  pressure, 

3  i  1  -I 

—  +  -  V<£  •  V<£-  rixr  •  Vtf>+  -(rixr)-(rixf) 

at  2  2 

+  -  +  gh+  (Oxd  +  nxrixd)  r  =  f{t)  {1.4) 
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However  reducing  the  problem  to  two  spatial  dimensions,  as  in  Figure  1.1,  where  the 
vessel  rotates  about  some  fixed  point  P,  results  in  the  set  of  equations  given  in  Section 
1. 


Appendix  II 
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In  this  appendix  the  coefficients  for  the  second  order  velocity  potential  and  wave 
height  for  the  elliptic  standing  wave  are  given. 
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Appendix  III 


In  this  appendix  a  theorem  for  the  solvability  of  the  inhomogeneous  wave  equa¬ 
tion  is  given.  A  variation  of  this  solvability  condition,  when  m=l  and  the  boundary 
conditions  are  dirichlet,  was  first  derived  by  Keller  and  Ting[8]  and  used  by  them  to 
prove  non-existence  of  periodic  solutions  for  a  class  of  nonlinear  wave  equations.  It  is 
used  here  to  prove  existence,  albeit  formal  existence,  for  the  nonlinear  wave  equation 
governing  standing  waves. 

Theorem:  Given  the  inhomogeneous  wave  equation, 


2d2T  d2r 

J°lh*  ~  fa* 


=  F(z,t),  0  <  x  <  1,  -oo  <  t  <  oo 


dr,  ,  dr,  . 

-(0,t)  = -(1,0  =  0 


(iii.i) 


(HI- 2) 


where  u>o  =  am  =  mn,  m  is  any  integer  greater  than  zero,  T z,  T t,  and  F  are  periodic 
in  time  with  period  2n,  F{x,t )  is  an  even  function  of  x,  and  in  addition  satisfies 


F(x  +  l,t  +  7r)  =  F(x,t), 


(III.  3) 


then  there  is  a  solution  to  the  inhomogeneous  problem  if  and  only  if 

J 0  I2am’  2  j 


(III.  4) 


Proof:  Every  solution  of  the  homogeneous  problem 
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which  is  27t  periodic  in  time  may  be  expressed  as 
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Multiply  (III.l)  by  ^(x,<)  and  integrate  from  0  to  1  in  x  and  0  to  2 it  in  t.  Integration 
by  parts  and  the  end  and  periodicity  conditions  result  in  the  left  hand  side  becoming 
zero,  then 

f  f  '&(x,t)F(x,t)dxdt  =  0  (III. 8) 

Jo  Jo 

For  the  higher  modes  the  integrand  may  be  expanded  in  x  taking  advantage  of  the  even 
extension  of  F(x,  t)  to  all  x, 

—  f  f  ty(x,t)F(x,t)dxdt  —  0  (///.9a) 

m  Jo  J - ~ 

when  m  is  even  and 

2  7T  i*  m  i, * 

-  I  I  *(x,t)F(x,t)dxdt  =  0  (III.9b) 

m  J0  /_  '-—JLL 

when  m  is  odd. 


Now  a  change  of  variable  is  made  to 


$  =  t  +  amx  (///.10a) 

X  =  t  -  ctmx  (III .106) 


The  region  of  integration  in  the  x  plane  for  the  first  and  second  modes  is  shown  in 
Figure  III.l.  The  higher  modes  have  regions  of  a  similar  nature.  With  the  change  of 
variables  (III.  10),  and  use  of  (III.3),  and  the  fact  that  the  integrand  is  similar  in  any 
period  parallelogram  in  Figure  III.l,  the  integral  (III.9)  becomes 


f '  f "  N + [fe:- : - 0  'mu> 

which  may  be  written  as  the  sum  of  two  integrals  f(£)F  and  f(x)F .  In  the  first  of  these 
the  variables  of  integration  are  interchanged.  Since  F  is  even  in  x  the  sum  becomes 


*L  mL  F {,um 

As  a  consequence  of  the  arbitrariness  of  f(x),  this  implies  (III.4)  and  proves  the  theo¬ 


rem. 


Figure  III.l  The  range  of  integration  in 
the  transformed  domain  for  mode  1  and 
2.  The  region  with  horizontal  stripes  is 
for  m=l  and  the  vertical  striped  region 
is  for  m=2. 


Appendix  IV 


In  this  appendix  the  details  of  the  Jacobian  Elliptic  functions  necessary  for  the 
construction  of  the  cnoidal  standing  wave  solutions  are  given.  They  are  taken  from  the 
Handbook  of  Elliptic  Integrals  [l]. 


The  complete  elliptic  integral  of  the  first  kind  is 


K(l 


y/l  -  x2sin2f 


(7K1) 


where  in  the  standing  wave  problem  K(f ;x)  =  n.  The  incomplete  elliptic  integral  of 
the  second  kind  is 


E(x;x)  =  J  |l  -  k2  + /c2cn2(f;/c)jdf 


(IV.  2) 


which  when  K(  >c)  =  ir  has  the  completed  form 


E(7r; /c)  =  J  yjl  -  k2  sin2  f  df 


(IV.  3) 


The  Jacobi  Zeta  function  is 


Z(x; k)  =  E(x;x)  - 

2 1  K) 


(/V.4) 


which  has  the  complete  Fourier  series 
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where  q  -  exp(-7rK(|; »c')/K(|; «))  and  k'  =  \/l  -  «2  and  the  cnoidal  function  has 
the  Fourier  series 
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cn2(x; x)  =  cos[j77^r] 
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for  n  >  0 
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The  constants  d j,  d2  used  in  the  expression  for  w2  in  Section  4  are 


d\  =  12/c2  — - —  45(1  -  k2)  +  8A(2k2  -  1) 
B 


( IV.8a ) 


d2  =  -6k2  ^  -  2(2k2  -  1) 
B 


(IV. 8b) 


where  A  and  B  are  defined  in  Section  4.  The  integrals  used  in  the  first  order  solution 
as  well  as  u>2  are  now  given  where  Cp  =  /02am  cnp(f;  /c)df  and  K(|;  k)  =  ir 
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